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Optomechanical devices in which a flexible SiN membrane is placed inside an optical cavity 
allow for very high finesse and mechanical quality factor in a single device. They also provide 
fundamentally new functionality: the cavity detuning can be a quadratic function of membrane 
position. This enables a measurement of "position squared" (x 2 ) and in principle a QND phonon 
number readout of the membrane. However, the readout achieved using a single transverse cavity 
mode is not sensitive enough to observe quantum jumps between phonon Fock states. 

Here we demonstrate an x 2 -sensitivity that is orders of magnitude stronger using two transverse 
cavity modes that are nearly degenerate. We derive a first-order perturbation theory to describe 
the interactions between nearly- degenerate cavity modes and achieve good agreement with our 
measurements using realistic parameters. We also demonstrate theoretically that the x 2 -coupling 
' ff l , should be easily tunable over a wide range. 

On 

i 1- I. INTRODUCTION 

In quantum mechanics a system's behavior is not independent of how it is measured. As a result, the readout used 
in an experiment must be tailored to the phenomena of interest. Likewise, for a given type of readout not all quantum 
C$ 1 effects are observable. 

Experiments on mechanical oscillators have to date used readouts that couple directly to the oscillator's displace- 
ment. The most common example is an optical interferometer in which the oscillator serves as one of the interfer- 
ometer's mirrors. In such a system the phase (j) of the light reflected from the interferometer is proportional to the 
mirror's displacement x. An oscillator that is subject to continuous monitoring of x is predicted to show a number of 
^ ' striking quantum features, including the standard quantum limit of displacement detection. 1 Additionally, the linear 
coupling between x and <p can be used both to laser-cool the oscillator (perhaps eventually to its ground state)^^ 
and to squeeze the light leaving the cavity^ The connection between the readout of the mechanical oscillator and 
its manipulation highlights the fact that these are two aspects of the same optomechanical coupling. 

In a recent paper 7 it was shown that a modest rearrangement of the usual optomechanical setup can realize a 
fundamentally different type of readout. When a nearly-transparent dielectric membrane is placed inside a cavity 
formed by two fixed, macroscopic mirrors, the phase of the light reflected from the cavity can be adjusted so that 
it is proportional either to x or to x 2 . The quadratic coupling occurs when the membrane is placed at a node (or 
anti-node) of the intracavity standing wave. In such a situation the membrane is at a minimum (maximum) of the 
optical intensity, and so detunes the cavity resonance by a small (large) amount. As the membrane moves in either 
direction it encounters an optical intensity that is larger (smaller) by an amount quadratic in its displacement (to 
lowest order), and hence detunes the cavity by an amount which is also quadratic (to lowest order) in the displacement. 
If on the other hand the membrane is originally placed at a point which is neither a node nor an antinode, the cavity 
detuning is (to lowest order) linear in the displacement. 

Mechanical oscillators coupled to an x 2 -readout have been discussed theoretically for some time. It has been shown 
that such a readout, coupled to a mechanical oscillator inside a sufficiently high- finesse optical cavity, can in principle 
provide a quantum nondemolition (QND) measurement of the energy (or equivalently the phonon number) of the 
mechanical oscillator. 8 With a sufficiently sensitive x 2 -readout it should be possible to observe, in real time, the 
individual quantum jumps of the mechanical oscillator. This is in contrast to an oscillator coupled to an x-readout, 
in which the repeated measurements of the oscillator's position extract information which prevents the oscillator from 
remaining in an energy eigenstate. This is because the quantity x does not commute with the oscillator's energy, 
whereas the quantity x 2 does (at least in the rotating-wave approximation, whose validity is ensured by the cavity's 
high finesse). 8 

Although the x 2 -readout demonstrated in reference [7] represented a major advance towards realizing the goal of 
QND measurements of a mechanical oscillator's energy, the strength of the x 2 -coupling was insufficient to realize such 
a measurement in practice. This is because for a low-reflectivity membrane the scale of the x 2 -coupling is ~ 1/A 2 , 
where A is the wavelength of the light (the cavity detuning oscillates each time the membrane is displaced by A/2). If 
the membrane's (field) reflectivity r approaches unity, the finesse of the "half-cavities" on either side of the membrane 
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FIG. 1: (a) Schematic of our apparatus: A flexible SiN membrane mounted on a motorized tilt stage at the center of a Fabry- 
Perot cavity is coupled to the cavity's optical modes via radiation pressure. Piezoelectric actuators between the mounting plate 
and membrane enable displacements along the x-axis. (b) Simplified diagram of the cavity and membrane. The cavity length 
is L — 6.7 cm and the end mirror radius of curvature is 5 cm. 

begins to increase, and the curvature of the cavity detuning (and hence the strength of the x 2 -coupling) increases, 
diverging for r — > l. 7 However the technical challenges involved in combining a high reflectivity mirror and a high- 
quality mechanical oscillator into a single element have proven considerable, so it would be highly advantageous to 
find a strong x 2 -coupling which does not require a high-reflectivity membrane. 

In this paper, we describe a new means for generating a strong x 2 -coupling in this type of device. We show that the 
optical cavity's full spectrum of transverse modes contains many near-degeneracies, and that near these points the 
cavity's resonance frequencies display an avoided-crossing behavior as a function of the membrane displacement. This 
leads to a detuning proportional to x 2 , but with a scale set by the symmetry-breaking aspects of the cavity /membrane 
geometry rather than the wavelength of light. We develop a perturbation theory that allows us to calculate the 
membrane-induced cavity detuning, and find that the x 2 -coupling at these avoided crossings can be made orders of 
magnitude stronger than realized in earlier work. We compare these calculations to measurements and find quantitative 
agreement, indicating that the single-phonon QND measurements proposed in [7] may be feasible even with a low-r 
membrane. 



II. OBSERVED EFFECT OF MEMBRANE ON EMPTY-CAVITY MODES 

Our experimental setup is shown in Fig. [1] and has also been described elsewhere 7,10 . A flexible silicon nitride 
membrane (1 mm x 1 mm x 50 nm thick) is situated near the waist of a high- finesse Fabry-Perot cavity so that its 
normal vector is roughly parallel to the cavity's long (x) axis. The membrane acts as the micromechanical resonator 
and its deflection is coupled to the cavity's optical modes via radiation pressure. The two macroscopic end mirrors 
are held fixed by an Invar cavity spacer. A motorized tilt stage holding the membrane is mounted to the spacer, and 
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FIG. 2: Transmission through the cavity as a function of laser detuning and membrane position. The dominant signal 
corresponds to the TEMo,o (singlet) cavity mode. Dashed lines show the approximate position of the unperturbed singlet 
modes. We have labeled positions corresponding to a node and antinode of the upper singlet mode's electric field. At these 
points the detuning is proportional to x 2 . (inset) An infrared camera image showing the transmitted beam profile. 



two piezoelectric actuators are used to displace the membrane along x. 

We can begin to characterize the optomechanical coupling in this system by measuring the transmission through 
the cavity as a function of membrane position and laser detuning, as shown in Fig. [2j Here the laser is aligned so that 
the dominant transmission peak corresponds to the TEMq,o (singlet) mode, as confirmed by a camera monitoring 
transmission (inset). As the membrane moves along the longitudinal (x) axis, it perturbs the cavity resonance 
frequencies to lower values, producing a detuning that varies roughly sinusoidally with position. 

When the membrane is located at an optical node, the perturbation is minimal, and the detuning is quadratic in 
position. As a result, light leaving the cavity contains only information about x 2 . As discussed elsewhere 7 -, this may 
enable QND phonon number readout using the TEMo,o mode alone. Since the membrane is a thin (50 nm) dielectric 
(nsiN ~ 2), it is a very poor reflector (|r| 2 = 0.13 where |r| 2 is the power reflectivity). As a result the curvature of the 
detuning is small and the x 2 -sensitivity is weak. Practical estimates predict that in order to observe a phonon Fock 
state before it decays, the membrane reflectivity would need to be substantially higher, ~ 0.998 7 . This may represent 
the most difficult of the technical challenges to observing real-time quantum jumps of the membrane's mechanical 
energy. 

A promising solution to this problem lies in the interactions between different transverse optical modes. We can 
couple to and identify many more of the cavity's transverse modes by intentionally misaligning the input laser, as 
shown in Fig. [3^a). We have identified all of the visible bands, such as the {TEMi, , TEM ,i} doublet, {TEM 2 , , 
TEMi 5 i, TEMo,2} triplet, and so on up to the 13-fold degenerate (tridectet) modes. 

The different transverse modes cross each other as a function of position at several places in Fig. [3l Figure [3] also 
shows a close-up of the crossing between the singlet and the triplet with the membrane's normal vector (b) aligned 
with x (c) tilted around the ?/-axis by 0.4 mrad. Tilting the membrane as in (c) lifts the degeneracy of the triplet 
in a predictable way: modes extended the furthest in the z direction shift the most. As is evident from Fig. Etc), 
the crossing points between the singlet and the two even triplets (TEM2,o and TEMo,2) are avoided, meaning that in 
addition to perturbing the individual modes, the membrane also couples them. 

Most importantly, the quadratic detuning at the avoided crossing turning points is very strong. In Fig. [3^c), 
the curvature is 50 times stronger than at the single- mode turning points. This is the same curvature a membrane 
reflectivity of 0.994 ± 0.001 would generate using a single mode 7 . This is already very close to the QND target and 
we have not yet attempted to optimize the system*^ 
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FIG. 3: (a) Transmission spectrum with the input laser misaligned, plotted on a log scale to enhance the faint features. The 
degeneracies of the different transverse modes are labeled and solid lines are drawn over the singlet mode for reference, (b) 
Close-up of the singlet-triplet crossing point for the membrane aligned with its normal vector parallel to the cavity axis (x). (c) 
Singlet-triplet crossing with the membrane tilted about the 2/- axis by 0.4 mrad. The strength of the curvature at the marked 
gap corresponds to an effective membrane reflectivity of 0.994 ± 0.001 power. 



III. MODEL 



While the membrane perturbs the empty cavity modes by up to a quarter of a free spectral range, its fractional 
effect relative to the laser frequency is minute (~ 10 -6 ). We can therefore view the membrane perturbatively and 
develop a first-order theory to model the system, as discussed in the next section. We then outline a method by which 
to solve this problem analytically when the membrane is positioned near the cavity waist, and finally compare our 
results with measurements. 



A. First-Order Degenerate Perturbation Theory 



We start with the time-independent free-space electromagnetic wave equation 

V 2 ^+^0 = O (1) 
c z 

where uo is the angular frequency and c is the speed of light. As is drawn in Fig. QJb), we define the origin to reside 
at the center of our cavity with the x-axis pointing toward one of the (spherical) end-mirrors. Under these boundary 
conditions, a convenient set of (Hermite- Gaussian) ortho normal solutions is given by 9 

, = H rn (V2y/w)H n (V2z/w) c _ (y 2 +z 2 )/w 2 
3 wV7rL2 m + n - 1 m\n\ 

x e i(rn+n+l)^ e -ik(y 2 +z 2 )/2R e -ikx-il7r/2 ^) 

Here H m is the m th Hermite polynomial, k = uo/c is the wavenumber, w(x) = \/2(x 2 + xj^j/kxR is the width of 
the cavity mode at x (where xr = 2.351 cm is the Raleigh range and wq = 89.2 fim is the waist for our geometry), 
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L = 6.7 cm is the cavity length, m and n are the transverse mode indices, I is the longitudinal mode index, ^f(x) — 
i&n~ l (w 2 k/2R) is the Guoy phase shift, and R(x) = (x 2 + xj^j/x is the wave fronts' radius of curvature. For a 
standing wave in our cavity, Re(^) is proportional to the electric field amplitude, and the In/ 2 term ensures that 
each longitudinal mode's electric field is zero at the end mirrors. The prefactors ensure that the inner product 
J dx J dy J dzRe(0i)Re(0j) = Sij. 

As shown in Fig. QJb), we represent the membrane as a block of refractive index nsiN ~ 2 and thickness t = 50 nm, 
centered at position xq. This modifies the speed of light in this short region, so that the wave equation in the cavity 
becomes 



uo 2 



V> + ^(1 + V(x - x c ))^j = (3) 

where V(x — x c ) = (n| iN — 1) (Q[x — (x c — t/2)] — Q[(x c + 1/2) — x]) and O the Heaviside step function. For an 
"aligned" membrane (i.e. flat in the y-z plane) x c = xq is constant. To incorporate tilt into the model, let x c — 
xo + ot y y + a z z where a y and a z are the small rotations about the z and y axes, respectively. 
The perturbed modes i/j can be expanded in terms of the empty-cavity modes: 

l/j = Ci0i + C 2 02 + C 3 (j) 3 + ... (4) 

where the c's are constants. We wish to study the region of near-degeneracy shown in Fig. Hfb-c), between the l-th 
longitudinal singlet mode {m = n = 0) and the three (I — l)-th triplet modes (m + n = 2), so we make the assumption 
that i/j is composed mostly of these four empty-cavity modes 

l/j = C S (j) S + Cy(j)y + C a (j) a + C Z (j) Z + ^2 € J ( t ) j ( 5 ) 

where the indices 5, y, a, and z refer to the singlet, the triplet widest in the y direction (m = 2, n = 0), the 
antisymmetric triplet (m = 1, n = 1), and the triplet widest along f (m = 0, n = 2), respectively. The last term is a 
summation over all remaining modes, and its contribution is assumed to remain small (ej <C 1). We also assume ip 
will have a new eigenvalue uo 2 /c 2 = n that is not very different (i.e. within a fraction of a free spectral range) from 
any of the unperturbed eigenvalues of the four contributing 0's. Substituting this into Eq. [3l 

(V 2 + (1 + V>)(C S (/> S + Cy<f>y + C a (t) a + C Z (j) Z + = - 

If we now take an inner product of this equation with each of the four empty-cavity modes and divide through by ft, 
we obtain four new equations 

(1 - Ki/K)Ci + V is c s + V^Cj, + Fi a c a + V iz c z + ^ Vi^-ej = (7) 

where the index i is s,?/, a, or z and V^j is the inner product of the z-th and j-th mode with V(x — x c ). The inner 
products Vij involve an integral over thickness t between two modes that are normalized over length L and are small 
(of order (n| iN — l)t/L ~ 2 x 10 -5 or less), so the last term in Eq. [71 can be ignored. We can further simplify by 
writing n = n s (l + 5) and K y , a ,z = + g) where 5 is the fractional change due to the membrane and g is the 
(constant) fractional separation of the unperturbed singlet and triplet bands due to the Guoy phase. Both S and g 
are of order 10 -5 . To first order, the remaining equation can be written as a matrix 
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Solving this eigenvalue problem for S in terms of g and the V's is straightforward and, though time-consuming, it is 
also easy to numerically compute Vij . The problem is in principle solved, and the result of such a calculation is shown 
in Fig. [U Computation time can also be reduced by assuming the membrane is an infinitesimally thin sheet (also 
plotted), but even a small finite thickness of t = 50 nm produces a noticeable effect. 



B. Analytical Solution Near the Cavity Waist 



We can also solve the inner products Vij analytically by using an approximate form for the unperturbed modes 
near the cavity waist. We do this by expanding 1/R(x), w(x), and *ff(x) in Eq. [2] in terms of the small parameter 
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FIG. 4: Comparison of numerical (dashed) and analytical (solid) results near the singlet-triplet crossing points. The red curves 
correspond to the thin-membrane (delta function) limit, and the black lines include membrane thickness. For this plot, xo = 
500 /iim, a z = 0.4 mrad, and a y = 0. 
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Substituting this into Eq. [2] yields 



H m (V2y/w)H n (V2z/w) 
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We can now use this to estimate the inner products Vij . We will not specify which modes are under consideration, so 
the results of this section may be applied to any set of nearly-degenerate cavity modes. 
Before we attempt to solve these integrals, first note that 
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Re^VRe^j) = Re 



(13) 



which simplifies the calculation. The first integral is by far the most challenging and so we outline its solution here. 

First, note that the singlet and triplet modes have slightly different unperturbed fc's and w's. If i = s and 
j G {y, a, z}, then kj = fc$(l + g) and wqj w wo^l — g/2). By defining A = 1 + g/2 ~ 1, we can easily keep track of 
this difference to first order (and A = 1 if the modes belong to the same degenerate manifold). If we plug Eq. [12] into 
the first integral of Eq. [131 make the substitutions y — > wyj\p2A^ z — > wz\p2A and x xrA, we have 



l e -Hh+h)^/2 



I dy I dz L 
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(14) 
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with 

_ (t/cosa)(n| iN - 1) 



, = (15) 

7rLy/2 mi+ni+m j +n jni\mi\nj\mj\ 

A t = t/xReosa (16) 

A c = x c /x R = A + l3yy + l3 z z (17) 

pfo,*) = H mi (y/VA)H mj (VAy)H ni (z/^ (18) 

if = 2^X^-^ + ^+771^+77,^+2). (19) 



Here A £ is the dimensionless membrane thickness corrected for tilt a = y cty + o^, and we have allowed the position of 

the membrane center x c to depend on y and z through small (rescaled) tilts in both directions, /3 ViZ = a y ^ z Wi/ xr\/2A. 

The membrane is 20 times thinner than the free-space wavelength, but it noticeably affects the cavity modes, as is 
evident in Fig. 0J We can approximate the integral over A by noting that for a smooth function /(#), 

x +5x/2 -i 

f(x)dx = Sxf(x ) + 5x 3 —f"(x ) + 0(Sx 5 ) (20) 



xq-Sx/2 24 



Applying this to our integral, 



HfaV « P Te-^ li+l ^ 2 JJ dy dz p(y,z)e- y e~ z 

x e -i[(A +[3 y y+[3 z z)(y 2 +z 2 )} 

x e -i[(A +/3 y2 y+/3^)K] (2\) 



with T = 1 — A 2 if 2 /24. It should be noted that K ~ 250,000 is very large, and so when estimating the thickness 
correction T from Eq. O we ignored several terms smaller than A 2 if 2 /24 ~ 10 -3 by a factor of K or more. 

The exponent in the second line of Eq. [21] contains only small quantities, so we can simplify this term by making 
the expansion e le « 1 + ie. Then we complete the square for y and z in the remaining exponential and make the 
variable change y — > y — i/3 y K/2 and 2 — > 2; — i/3 z K/2. If we then define the (analytically soluble) integral 

= Jdx(x- if3K/2)"e- x2 H ni ( - ~^ /2 )H nt ((x - i0K/2)VA) (22) 
and the shorthand T qp = ^ x m.^nf , Eq. |2"T1 becomes 



^ /// ~ Pore- iJrA °- i ('« +, i> ,r / 2 e- 



x [r 00 - % ((r 20 + r 02 )A + (r 30 + r 12 )^ + (r 03 + r 21 )&)] . (23) 

Applying a similar method to the second half of Eq. and with definitions K' = mj — mi + rij — rti — ghxR (much 
smaller than K) and T' qp = £^f.^ £nfnf~ , it can be shown that 
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jV « p oe -iK ^ a -i{U-l j)-K/2 

x [4 -^'(4A + r' 10 /3 y + r^)]. (24) 



Equations [23j [24] and [13] represent a very accurate analytical approximation of the inner products Vij for small 
displacements (relative to xr) from the cavity waist. These results are also plotted (with and without the thickness 
correction) in Fig. [4] of the previous section as solid lines. In practice, the agreement for our setup is excellent as long 
as \x c \ < 1 mm. As expected, the approximation is nearly perfect at the waist and breaks down as xq approaches xr. 
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C. Discussion and Comparison with Data 

We can gain insight into our system from the analytical results. By ignoring the off-diagonal terms in Eq. [8] (i.e. 
ignoring avoided crossings), the eigenvalues simplify substantially, and the detuning is given by 



Au;, 



*(rcjiN ~ 1) (1 



uq 2L 
coo 2L 



Tcos((2to-1)A )) (25) 



(1 + T cos((2k y x R - 3) A )) + g (26) 



for the singlet and triplet modes respectively. These equations represent two sinusoidal bands oscillating (with opposite 
sign) below their unperturbed detunings, with peak-to-peak amplitudes T(n| iN — l)t/L ~ 27% of the free spectral 
range, and separated from each other by the appropriate Guoy spacing. Applying this method to the other transverse 
modes, we can generate the entire band structure shown in Fig. [3](a), and for nsiN = 2.0 the agreement is essentially 
perfect. 

The spatial period of the singlet band is slightly smaller than the triplet, and so it should oscillate a little faster as 
a function of A . Though subtle, we do observe a phase difference between the bands (this is somewhat more visible 
in Fig. [3] when comparing the singlet and the nonet band), and we can use this phase difference to estimate the 
membrane's displacement xq. In the time it takes to raster the band structure in Fig [3j drifts in the piezos and laser 
frequency cause noticeable distortions of the bands on the scale of this phase difference, but nonetheless by fitting the 
neighboring singlet and triplet sinusoids from Fig. [3^a) we estimate a positive (as defined by the x-axis) displacement 
of 300 ± 100 /im from the waist. 

It is also worth noting that the effect of finite membrane thickness is to wash out the oscillations (i.e. T decreases 
from unity as t increases). This makes sense qualitatively because when the membrane is positioned at a node, the 
electric field is not zero everywhere inside it and so there will still be a small negative perturbation at the top of 
the band. Similarly, when the membrane is at an antinode, the field is not maximal everywhere inside it and so the 
perturbation is not as strong. Optical losses inside the membrane mean that even when positioned at a node the 
finite thickness will put an upper bound on the finesse this system can achieve. We have, however, already observed 
a finesse of 150,000 with the membrane inside the cavity^. 

If we now look at the mode-coupling terms V^-, we can gain some insight into the avoided crossings of Fig. [3jb-c). 
First, all of the off-diagonal terms involving the antisymmetric mode (TEM^i) are identically zero (even with 
in this approximation. This is a reflection of the fact that the TEM^i mode is an odd function in both the y and z 
directions, while the other three modes are even. The integrals across the membrane therefore all involve a function 
that is approximately odd and vanish. Hence there should be no avoided crossing between TEMo,o and TEM^i to 
first order, which agrees with all of our observations (see Fig. Etc), for example). 

The other off-diagonal terms are not zero (thankfully), and the result is again relatively simple if we keep the 
membrane aligned (i.e. a = 0). 



V sy = V sz « -A * (ns ' 2 N £ 1) '/ ..,> [((k s + k z )x R - 4)A ] . (27) 



When the membrane is aligned, the interaction is proportional to Ao times a term that oscillates with a period close to 
that of the bands (though since the bands always cross each other at roughly the same phase, this term will modulate 
the coupling slowly as a function of Ao). Following this backward through the calculation in the previous section, we 
see that it arises from our expansion of the finite radius of curvature R. So (perhaps not surprisingly) the interaction 
between these modes arises from the mismatch between the curved wavefronts and the fiat membrane^ 

This is an encouraging result because it implies a strong degree of tunability in the avoided gap and hence the 
-sensitivity. Figure [SJa) shows our calculation of the gaps at each of the four avoided crossings as a function of 
displacement out to 500 //m from the waist. In this plot we show results for both the aligned case (dashed lines) and 
for a tilt of 0.25 mrad (solid lines). For the aligned case, the gaps collapse onto two similar curves (as they must 
by symmetry), and when the membrane is tilted, the top two gaps (TR and TL of the inset) move to larger initial 
values. 13 The top right (TR) gap is a more interesting function of position, as it is tunable through zero at finite 
offset. If instead we fix the position of the membrane, we should also be able to tune the TR gap over a wide range 
(including zero) with tilt, as shown in Fig. [5fb). Note the large quantitative and qualitative differences between the 
TR and TL gaps can help calibrate the magnitude as well as the direction of the membrane displacement if it is not 
already known. 

As mentioned in section[IIl vibrations in this apparatus preclude reliable determinations of gap size and x 2 -sensitivity 
for very small gaps. Nonetheless we have observed some smaller gaps as we tune the membrane tilt, two of which are 
shown in Fig. [6j In this data set, displacement noise due to ambient vibrations, coulomb forces on the membrane 
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FIG. 5: (a) Dependence of the four avoided crossing gaps on membrane position, fixing tilt at (dashed lines) and 0.25 mrad. 
(inset) Plot of the mode detuning versus position near the singlet-triplet crossings. The four gaps in (a, inset) are labeled for 
reference, (b) Dependence of the four gaps on membrane tilt, fixing the position at 300 and 800 /mi. The TR gap should be 
adjustable over a wide range. 



and/or piezo noise is quite evident. Further, during the left-to-right rastering of these data sets (acquired over ~ 10 — 20 
minutes each) the laser temperature varied enough to cause a systematic detuning and sheer the data vertically. It is 
a large effect in this data set, and we are studying ways to compensate for it. 

In a given frequency sweep (vertical trace), however, the time it takes to traverse one of these gaps is roughly 
a millisecond; such a single-shot measurement of the mode spacing should therefore be much less susceptible to 
vibrations and drift. If we therefore record the smallest spacing in Fig. [6l we can put a lower bound on the detuning 
curvature using the form detuning takes near an avoided crossing, \J (ax) 2 + (A//2) 2 where a is the asymptotic slope 
and Af is the gap (both of which we estimate from Fig. [6]). Doing so yields a lower bound on the effective membrane 
reflectivity of |r| 2 > 0.992 ± 0.004 for the TR crossing and |r| 2 > 0.9989 ± 0.0005% for the BR crossing (the sharper 
curvature of BR reflects the larger asymptotic slope a). This estimate is still subject to vibrations above a few 
kilohertz, which we have not characterized. On the other hand, when we fit the curvature explicitly as in section [Til 
even lower- frequency vibrations (i.e. anything above about 0.1 Hz) can wash out sharp curvature, so that technique 
represents a very conservative estimate. 

Figure [5] implies we can use the TR and BR gaps to estimate the position of our membrane relative to the waist. 
Figure shows the aligned singlet-triplet crossing data from Fig. [3f b) along with curves generated by this model 
(for a = 0) at several different membrane positions. We estimate the membrane's displacement from the waist to 
be about 550 jam here, which is a reasonable value for our apparatus and is in rough agreement with our previous 
estimate based on the horizontal offsets in the various bands shown in Fig. [3l 

We can further check the model for consistency by studying the interplay between tilt and the lifting of triplet 
degeneracy far from a crossing. It is relatively straightforward to show that this scales as a 2 for small a. The triplet 
splitting should also be quite insensitive to membrane position so we can use it to estimate the membrane's true tilt or 
even align the membrane^ Figure [8] shows the data from Fig. [3](c) along with the analytical result for xq = 325 /im 
and a z = 0.395 mrad. We obtain these parameters by first adjusting a z until the triplet splitting is correct and then 
varying xq to match the avoided crossings. We have also plotted the result for displacement in the opposite direction, 
which essentially amounts to comparing our data with BL and TL in Fig. [5j The fit does not agree with the data here 
or at any other negative value of xq. The smaller TR gap therefore confirms the sign of our membrane displacement 
(and it is the primary reason we chose to study the right side). 

We have performed similar analysis at several different values of the tilt stage's motor position, and these are 

summarized in Fig. [9l Assuming there is a small constant tilt a y , we can fit this data with the form a = (aq z ) 2 + a 2 

where a is a mechanical conversion factor between motor position q z and tilt. From the fit a y = 0.16±0.01 mrad and 
a = 0.0756 ± 0.0001 mrad//im. From the length of the tilt stage lever arm (12.7 mm) alone we estimate a = 0.0787 
mrad//im, implying a calibration error of ~ 4%. We have plotted the expected result for the same a y using a 12. 7- mm 
lever arm for reference. 

It is also important to note here that in the model we reproduce the ordering of the triplet modes: as we rotate the 
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1.562 1.568 1.574 1.580 1.586 1.592 

Ax IK 



FIG. 6: Transmission data at large membrane tilt, a z ~ 0.65 mrad. The dashed line corresponds to the data shown in the 
inset, (inset) Single-shot measurement of transmission versus laser detuning. 




1.574 1.577 1.580 1.583 1.586 1.589 



Ax /A 



FIG. 7: Analytical model plotted on top of transmission data for the aligned membrane. Here we show the analytical results 
for the membrane situated at 200, 550, and 1000 /mi from the cavity waist. 



membrane about the y axis, the modes most extended in the z direction move the furthest. 

The model is in reasonable agreement with the data thus far, and it implies that in future experiments we should 
be able to tune the x 2 -sensitivity to essentially any desired value. This could be a very important tool in our attempt 
to perform QND measurements of a single phonon. 



IV. SUMMARY/OUTLOOK 

In this paper we have demonstrated that a SiN membrane can couple two nearly-degenerate transverse optical cavity 
modes, generating an avoided crossing and a cavity detuning that is strongly quadratic in membrane displacement 
x. Without optimizing the system, we have shown that this x 2 -dependence (which is tunable over a wide range 
via membrane tilt) can be as strong as that generated using a single cavity mode and a membrane of reflectivity 
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1.568 1.574 1.580 1.586 

Ax IK 

FIG. 8: Model plotted on top of the transmission data for the tilted membrane. Here we show the analytical results for the 
membrane situated at +325 and -325 /mi from the cavity waist. 

|r| 2 > 0.9989 ± 0.0005. This means it might still be possible to perform QND measurements of phonon number 
in a membrane of modest reflectivity (i.e. |r| 2 ~ 0.13). We also derived a perturbative model of the system that 
quantitatively agrees with observations and further predicts the x 2 -strength should be tunable to arbitrary strength 
through mm-scale membrane displacements. 

These results should be taken with the caveat that the sharp avoided crossings described above occur when the 
membrane is not at a node of the intracavity field. As discussed previously^^, this means that the optical loss in the 
membrane will limit the maximum cavity finesse. Whether or not the effect of this reduced finesse can be offset by the 
very strong quadratic coupling or reduced optical loss (e.g., via improved membrane materials or further engineering 
of the cavity modes) remains to be seen. 
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FIG. 9: Plot of tilt estimated from triplet splitting versus our tilt stage motor position. The red curve is a fit allowing the 
motor's linearity constant and misalignment a y to float and the blue curve is the expectation for the same a y and determining 
the linearity constant from the tilt stage geometry. 



